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We study hydrogen doping effects in an iron-based superconductor LaFeAsOi_ y by using the 
first-principles calculation and explore the reason why the superconducting transition temperature 
is remarkably enhanced by the hydrogen doping. The present calculations reveal that a hydrogen 
cation stably locating close to an iron atom attracts a negatively- charged FeAs layer and results in 
structural distortion favorable for further high temperature transition. In fact, the lattice constant 
a averaged over the employed supercell shrinks and then the averaged As-Fe-As angle a approaches 
109.74° with increasing the hydrogen doping amount. Moreover, the calculations clarify electron 
doping effects of the solute hydrogen and resultant Fermi- level shift. These insights are useful for 
design of high transition-temperature iron-based superconductors. 

PACS numbers: 74.25.Jb, 74.70.-b, 71.15.Mb 



Since the discovery of the high-temperature supercon- 
ductivity in LaFeAsOi-^F^ [l|, iron-based superconduc- 
tors have drawn wide attention due to not only their high- 
temperature transition but also peculiar superconduct- 
ing properties. According to the accumulated insights, 
most of iron-based superconductors require chemical dop- 
ing like substitution as LaFeAsOi-^F^ or deficiency as 
LaFeAsOi-^ for their optimal high-T c superconductivity 
0| . In addition, most of their mother compounds ex- 
hibit the stripe- type antiferromagnetic ordering just after 
the orthorhombic transition on cooling the temperature. 
There are several similarities with cuprate high-T c super- 
conductors but some original particularities. Therefore, 
more systematic research is now still in great demand. 

Recently, Miyazawa et al.,[3] reported a fascinating 
fact that the transition temperature T c of LaFeAsOi-^ 
is enhanced from 28 K to 38 K by doping the hydro- 
gen in its high-pressure synthesis. It surprisingly cor- 
responds to more than 30% of T c enhancement, which 
is unlikely in other superconductors. This new synthe- 
sis scheme is quite different from the previous chemical 
doping such as substitutions and deficiencies, because 
the enhanced T c is much beyond the optimum reach of 
so-called La-1111 compounds. Therefore, the hydrogen 
doping effect should be an important key to the super- 
conducting mechanism and T c enhancement in iron-based 
superconductors. In this paper, we study the hydrogen 
doping effects through systematic first-principles calcula- 
tions. The highlights of this paper are to clarify at which 
the doped hydrogen atom prefers to stay inside the crys- 
tal unit and how the doped hydrogen atoms deform the 
crystalline structure and resultant electronic one. Most 
of the calculated results are consistent with measurement 
data, while they contain novel insights never experimen- 
tally accessible. 

Before starting the calculation study of the hydrogen 



doping effects, there are a few questions to narrow the 
point. The first one is whether the hydrogen atoms re- 
ally keep their positions inside the crystal. The nuclear 
magnetic resonance measurement on the hydrogen (H- 
NMR) confirmed the existence of doped hydrogen atoms 
d, [H[ not at any surfaces and boundaries but inside the 
crystal. This indicates that there are stable locations in- 
side the crystal for the hydrogen atom. The next one 
is internal locations of the doped hydrogen atoms. One 
naively expects that doped H atoms fills O vacancies. 
Contrary to this conjecture, Fe-NMR verified that the 
hydrogen atoms partly locate close to Fe atoms [5[ . How- 
ever, their exact locations have not been known yet be- 
cause of its technical difficulties. Then, we explore the 
exact position at which the hydrogen atom stably locates 
around Fe atom by first-principles calculations. Once one 
solves the above problem, the final goal is to clarify how 
the hydrogen atom raises T c . In this paper, we reveal 
some important roles of the doped hydrogen atoms. The 
first one is that the hydrogen atoms tend to shrink lat- 
tice constants like the replacement of Ln element in LnO 
blocking layer of 1111 compounds. The superconduct- 
ing transition temperature is known to be quite sensitive 
to the variation of the lattice constants. Moreover, it is 
well-known that the As-Fe-As angle, a and the height of 
an As atom from an Fe-plane, Has are closely correlated 
with T c @, 0]. In fact, we find that the doped hydrogen 
distorts the crystal structure more suitable for high-T c 
superconductivity. The second one is whether the hydro- 
gen doping supplies additional carriers. The calculations 
reveal that the hydrogen provides electron and the Fermi 
surfaces change their shapes in a large doping level. 

Let us briefly show the present calculation techniques. 
We examine hydrogen-doping effects in iron-based su- 
perconductors based on the density functional theory 
(DFT). The calculation package employed throughout 
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this paper is VASP Q, which adopts PAW method [To| 
and GGA [9] exchange-correlation energy. In order to 
examine how the doped hydrogen deforms the crystal 
structure, we prepare 2x2x1 supercell of (LaFeAsO)s as 
shown in Fig. [TJa). If one puts a hydrogen atom into the 
supercell, the chemical formula is given by LaFeAsOH^, 
in which x = 0.125. Accordingly, x — 0.25 when insert- 
ing two H atoms. Thus, we can control the doping level 
x, although the variation range is quantized. In order 
to find relaxed optimal positions for the hydrogen atoms 
inside the supercell, we optimize the crystal structure of 
LaFeAsOH^, so that the forces on all atoms become less 
than 0.01 eV/A. The parameters set in the DFT calcula- 
tions are given as follows, /c-points are taken 4x4x4, and 
the self-consistent loops are repeated until the energy de- 
viation becomes less than 10 -5 eV with the cut-off energy 
being 500 eV. 

Now, let us present the calculation results. First, 
we explore stable positions for a hydrogen atom inside 
LaFeAsO crystal unit. We initially prepare a few tens 
configurations, in which an H atom is randomly dis- 
tributed inside the LaFeAsO unit cell, and optimize their 
structures. Then, the doped hydrogen atom finds out a 
stable location. As a result, we find two typical locations 
for an additive H atom. The first type is around an alter- 
nately vacant As position which is mirror symmetric with 
a non- vacant As across the Fe-plane. According to the 
calculation result, the exact stable location slightly shifts 
from the vacant As position as shown in Fig. QJb). The 
second one is just below a La atom alternately symmet- 
ric across the O plane as shown in Fig. Die). The former 
case is found to be more stable by comparing their to- 
tal energies. This result actually agrees with Fe-NMR 
measurements [5j. However, the calculated energy dif- 
ference is 50 meV per LaFeAsOHo.125 and the hydrogen 
atom may partially occupy the second position as a quasi- 
stable location. 

Here, we turn to O vacancies. As mentioned above, H 
atoms are expected to enter O vacancies. In order to ex- 
amine this, we compare the total energies of the following 
two types of (LaFeAs^OyH structures. A configuration 
is that a H atom fills the O vacancy, and the other one is 
that it locates near an Fe atom without filling the O va- 
cancy. As a result, more stable configuration is found to 
be the former one. Therefore, we conclude that H atoms 
first tend to fill O vacancies. In this case, the H doping 
corresponds to just a substitution like LaFeAsO i-^F^, 
and it is found to be not relevant to more than 30 % 
T c enhancement. In fact, the Fe-NMR result [5] insists 
that H atoms exist near Fe atoms. This indicates that 
an excess portion of the doped H atoms located in FeAs 
layers. Hereafter, we focus on the H locations near Fe 
atoms as a new mechanism of T c enhancement. 

Let us concentrate on the location inside FeAs layers. 
As shown in Fig. [ljb), we find out that the most stable 
position shifts to (0.07,0.07) in the inner atomic coordi- 




FIG. 1: (a) 2 x 2 x 1 supercell of LaFeAsO. (b) Most stable 
structure of (LaFeAsO)sH, in which a H atom exists in the 
FeAs layer, (c) Metastable structure of (LaFeAsO) 8 H which 
contains a H atom around the LaO layer. 
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FIG. 2: Total energy per LaFeAsO Ho- 125 as a function of 
inner atomic coordinates of the H atom in (LaFeAsO) 8- Inset 
shows the position of the H atom. Inner atomic coordinate 
zh, which is along c-axis, is optimized. 



nate from (0,0) corresponding to just above an Fe atom. 
We carefully trace a variation of the total energy with re- 
spect to the position shift as (xh,^h,^h), in which once 
xh is fixed is optimized. Figure [2] shows the total en- 
ergy variation as a function of xu. At xu = 0, the total 
energy counterintuitively shows maximum. The present 
calculation reveals the reason as follows. The position 
xyl = is actually unstable due to the repulsion from 
the nearest La 3+ cation for the H + one. Then, it shifts 
to avoid the so-called Coulomb-repulsion pressure from 
the La 3+ cation. This result obviously indicates that 
Ln cation in 1111 crystals has an important influence on 
FeAs layer. The Ln cation and FeAs layer keep exquisite 
balance. On the other hand, we notice in Fig. [2] that the 
energy difference in the energy valley is not so large that 
the hydrogen atom can largely vibrate around the stable 
position. Such a localized dynamical effect is beyond the 
present study. 

Next, we examine the hydrogen-doping amount depen- 
dence. Since we now know the most stable location inside 



3 



500 
„ 400 
| 300 
% 200 
100 

( 

4.03 



LU 



°< 



°< 



(a) 



4.02 

4.01 

4.00 

3.99 
9.20 



J p5T 



1 I ^ 

1 • • • 




0.4 0.6 
Doping x 



CD 



D) 
C 
CC 
CO 
< 
I 

CD 
LL 
I 

CO 

< 




0.4 0.6 
Doping x 



FIG. 4: As-Fe-As angles as a function of doping x. Small 
squares and circles stand for a and f5 angles, respectively, 
which are depicted in the inset. Large squares and circles are 
the average of a and /3, respectively. Horizontal line denotes 
the regular tetrahedron angle, 109.47°. 



FIG. 3: Ef(x) (see Eq. ([T])) (a), lattice constants a (b) and c 
(c) as a function of doping x. 



the crystal unit, as initial configurations we put H atoms 
on the positions as shown in FigOJb). There are 1 to 
8 As atoms inside the present super-cell. This indicates 
that there are also 8 possible initial locations for an H 
atom. Then, the number of possible configurations are 
totally 2 8 = 256. Among them, several configurations 
are symmetrically equivalent to others, and the initial 
sets can be then reduced to 22 patterns, among which 
the number of doped hydrogen atoms is ranged from 1 to 
8. We calculate stable structures and total energies on 
these all 22 patterns. Figure [3^a) shows an energy Ef(x) 
as a function of the doping level x, which is obtained by 
dividing the number of H atoms in the cell by 8. The 
energy Ef(x) is defined as 



£(LaFeAsOEy -^(LaFeAsO) - -£(H 2 ), (1) 



where £■(•••) stands for the total energy of the corre- 
spondent chemical formula. As seen in Fig. [3ja), Ef(x) 
is positive and increases with x. One finds that all the 
present possible doping is energetically unstable com- 
pared to non-doped LaFeAsO. This is related to a fact 
that high pressure synthesis is needed for doping. The 
doping x variation of the lattice constants a and c are 
shown in Figs. [3jb) and (c), respectively. From both 
the figures, one finds that a decreases with increasing 
x. The observed a of H-doped LaFeAsOi-^ is 3.9938 
A[3[, which is distinctly smaller than that of the non- 
doped LaFeAsO X - y (4.0257 A) [4]. This fact is consistent 
with our calculations. On the other hand, the observed 
c counterintuitively decreases from 8.7190 A [4] to 8.6898 
A[3[ by H doping, although the lattice constant c in- 



creases with x in the present calculation. The present 
tendency contradicts the measurement data. We suggest 
a reason on the opposite tendency as follows. It is well- 
known that the hydrogen can easily diffuse in various 
compounds because the energy barrier for the diffusion 
is not so high as seen in Fig|2] This means that any hy- 
drogens do not always occupy the energy minimum point. 
For instance, if a hydrogen atom partially occupies the 
second energy minimum location as shown in FigH] (c), 
then c-axis constant may rather decreases, because the 
neighbor is negatively-charged O 2- . Although such sta- 
tistical analyses on various configurations including the 
diffusion effects are beyond the present work, the c-axis 
shrinkage clearly indicates that doped hydrogen atoms 
partly stay at the positions close to the anions. 

Figure |4] is the highlight of the present paper. Lee et 
al., |6[ pointed out that T c is closely correlated to the 
As-Fe-As angle, a (see the inset of Fig. HJ where f3 is 
noted to be also another equivalent index). T c takes the 
maximum value when a(/3) = 109.47°, at which FeAs4 
tetrahedron becomes a perfectly regular one. As seen in 
Fig. HI the averaged a(/3) clearly approaches to the op- 
timal value with increasing the hydrogen doping ratio x. 
The reason is given as follows. The hydrogen ion (cation) 
attracts negatively-charged FeAs layer, and then the in- 
plane lattice constants decrease. Consequently, the As 
height increases and the angle a(/3) decreases (increases). 

The next issue is an effect of the hydrogen doping on 
the electronic band structure. We display band disper- 
sions of H-doped LaFeAsO in Fig. [U It is noted that the 
current Brillouin zone (BZ) is folded to a quarter due to 
the use of 2 x 2 supercell. In the conventional unit cell, the 
iron-based superconductors generally have hole pockets 
around the T point (the center of BZ) and electron pock- 
ets around M point (the corner of BZ). In the present 
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FIG. 5: Band structure of LaFeAsOH^ of (a) x = 0.125 (b) 
x — 0.25. The Fermi energy is set to zero. 



zone, M point moves to T one by the zone folding, and 
both the pockets almost overlap around T one. Actually, 
both the pockets crosses in the case at x = 0.125 as dis- 
played in Fig. [5ja), i.e., the Fermi surfaces between the 
hole and electron pockets are well nested. Meanwhile, 
the hole pocket almost disappears at x = 0.25 as shown 
in Fig. [5fb). The reason is that the doped hydrogen 
atom becomes a cation and supplies an electron. Con- 
sequently, the Fermi level lifts up and the hole pocket 
disappears beneath the Fermi level. This indicates that 
too much hydrogen doping drastically changes the elec- 
tronic structure around the Fermi level and eventually 
diminishes the high-temperature superconductivity. 

Finally, we discuss the hydrogen doping effects to- 
gether with various knowledge on hydrogen additives in 
solid states. First, we stress that the doped hydrogen 
shrinks all the lattice constants although extra element 
additives normally expand lattice constants. One of the 
origin is that the hydrogen cation locating stably close to 
Fe atom attracts negatively-charged FeAs layer. Thus, 
the mother compound LaFeAsO crystal is found to be 
rather loosely packed in contrast to very low solubility 
of the hydrogen in pure Fe crystal. This fact is deeply 
relevant to the large pressure and substitution effects in 
1111 compounds. From this discussion, it is expected 
that the hydrogen doping effect on Sm-, Gd- and Nd- 
1111 compounds leads to further T c enhancement. Un- 
fortunately, such enhancement was not reported, e.g., in 
H-doped SmFeAsOi- 



However, in their com- 
pounds, it is noteworthy that the hydrogen is also soluble, 
because their FeAs structures and associated electronic 
ones are still controllable. Second, we point out that 
the doped hydrogen supplies an electron. Namely, the 
hydrogen doping lifts up the Fermi level and the excess 
doping loses the nesting situation which is regarded as an 
electronic-structure stage for high-T c superconductivity. 
This means that there is the best hydrogen doping ra- 



tio while further T c enhancement can be expected if one 
controls the carrier number by different ways from the 
hydrogen doping, because the angle a monotonically ap- 
proaches the best angle with the doping amount. Third, 
we would like to mention that the local distortion gener- 
ated by the doped hydrogen considerably varies inside 
the supercell unit and just the averaged distortion is 
preferable for the transition temperature enhancement 
as shown in Fig. |4j Accordingly, the emergent su- 
perconductivity seems to reflect not local atomic-scale 
but widely averaged distortions and resultant electronic- 
structure changes. Since high diffusivity of the hydrogen 
is largely suppressed in a low temperature range as the 
superconducting transition, it is a big mystery that the 
interstitial hydrogen close to Fe seems not to work as a 
pair breaking potential. 

We have investigated the optimized crystal structure 
of the hydrogen doped LaFeAsOH^ by using the den- 
sity functional theory. The calculation results indicate 
that the hydrogen atom prefers to stay close to Fe atom. 
In the case, the supercell calculations revealed that the 
averaged lattice constant a shrinks with increasing the 
doping amount, which is consistent with the experimen- 
tal results. Meanwhile, the lattice constant c expands 
within the present calculations, although it also shrinks 
in experiments. Therefore, it is considered that the hy- 
drogen atoms locates at not only the most stable point 
but also other minimum ones. The averaged structural 
distortions coincide with the correlation tendency of T c 
on the crystal structure. As a typical example, we ob- 
served that the As-Fe-As angle a approaches closer to 
109.47° by increasing the amount of the doped hydrogen. 
However, the doping shifts up the Fermi level, because 
the hydrogen also has the carrier doping effect. In con- 
clusion, the present first-principles studies clarified some 
fundamental aspects in the hydrogen doping beyond the 
findings obtained by experiments. We believe that their 
insights are essential to materials design for further high- 
Tc superconductivity. 
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